Pure Lattice Gauge Theory in the Expanding Universe 



(N 

o 

(N 
=3 



(N 



43 



> 
O 

On 1 

OS' 

od 

o 

(N 



Marton Trencscnf] 
Dept. of Physics of Complex Systems, Eotvos University, 
Pdzmdny Peter setdny 1/A, 1117 Budapest, Hungary 
(Dated: August 27, 2012) 

Using Monte Carlo methods, I study the thermodynamic properties of Z2 Abelian lattice gauge 
theory in flat, homogeneous, isotropic, expanding spacetimes characterized by the scale factor a(t). 
The presence of the scale factor introduces a spacelike domain wall inside the four dimensional 
spacetime lattice where the phase transition associated with deconfinement occurs. The resulting 
theory appears identical to introducing a time-dependent effective coupling j3 e g, and could serve as 
a simple model for the electroweak phase transition of the early Universe. 
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Lattice gauge theory, the study of gauge theories on 
discretized spacetimes is an important tool for the study 
of quantum theories such as QCD and the Standard 
Model. On the lattice path integrals associated with the 
theory become finite dimensional, making it is possible 
to approximate them on the computer using Monte Carlo 
methods. Lattice gauge theories were introduced by Wil- 
son's 1974 definition of a discretized action with the ap- 
propriate continuum limit Q. Creutz' series of papers 
in the following years investigated the thermodynamic 
properties of Zn, U(l) and SU(N) pure lattice gauge 
theories using computer simulations, and found that the 
theories contain phase transitions associated with decon- 
finement 0, ■ This demonstrated that gauge theories 
are a viable framework for the experimentally observed 
phenomenon of quark confinement. 

These, and subsequent investigations that followed 
were carried out in Euclidean spacetime. In the context 
of this letter, it is instructive to think of Euclidean space- 
time as a Wick rotated Minkowski spacetime, where Wick 
rotation is defined by the replacement t — > it. Wick rota- 
tion is necessary to get rid of the imaginary exponent in 
the original path integral, to arrive at a theory which can 
be simulated effectively on a computer. But the results 
of the computational experiments describe the physics of 
gauge theories in Minkowski spacetime. 

Standard cosmological theory, confirmed recently the 
Wilkinson Microwave Anisotropy Probe [(§ and the 
Sloan Digital Sky Survey Q, tells us that the Uni- 
verse is expanding, not static. On cosmological scales 
the Universe is described by the spatially flat, homo- 
geneous and isotropic Friedmann-Lemaitre-Robertson- 
Walker (FLRW) metric, which contains a scale factor 
a(t) to describe the expansion of space with cosmolog- 
ical time. 

In this letter, I define pure lattice gauge theories on 
Euclidean FLRW metrics and present results of compu- 
tational simulations using the simplest Z2 Abelian model. 
These demonstrate that the presence of the scale factor 
introduces a spacelike domain wall inside the four dimen- 



sional spacetime lattice where the phase transition asso- 
ciated with deconfinement occurs. This indicates that 
even simple lattice gauge theory models extended with 
the scale factor a(t) can qualitatively reproduce the elec- 
troweak phase transition of the early Universe. 

I formulate the theory on a four-dimensional hyper- 
cubic lattice. Associated with each link joining a pair 
of nearest-neighbor sites i and j is an element Uij of 
the group Z 2 = { + 1,-1}, oriented such that Uij = 
Ujl . The discretized action describing these spins is 
S = J2a where the sum is over all elementary space- 
time plaquettes, and S\j is the per plaquette action de- 
fined as 



S 



□ 



g\ QmiQw (1 - UijUjkUkiUu) 



for a plaquette in oriented in the \xv plane, where g a p 
is the Euclidean FLRW metric (500 = 1, 9u = a 2 (t)) for 
the chosen cosmology, and g = det(g Q ^) = a 6 (t). The 
extra g terms in this expression arc there to reproduce 
the Einstcin-Hilbert action in the continuum limit, which 
for flat spacetime with matter is Seh = / £■ \f\g\ d A x. 

To see this, note that the Lagrangian density for the 
Abelian gauge field = d^A^ 



d v A^ is 



C = --FuvF^ = --g ai g? s 



where the group elements are connected to the vector 



potential by E/j. 



JgoA f 



d where go is the coupling 



constant of the theory, /i is the direction of the bond 
specified by ij, and d is the lattice spacing. 

We insert the discretized action into a path integral to 
define a partition function at inverse temperature (3 = 
V9l 



Z 



{Uij} 



where the sum runs over all possible values of the link 
variables Uij. The order parameter used to search for 
phase transitions is the average action per plaquette 

E = (1 - UijUjkUkiUu) 
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FIG. 1. Results of Monte Carlo runs showing the order param- 
eter E for an effective 20 4 lattice with different cosmological 
scale factors. The black line corresponds to /3 c ff(i) = fie- The 
plots demonstrate that a cosmological phase transition occurs 
when this effective critical inverse temperature is reached. 



Unlike in traditional pure lattice gauge simulations where 
the order parameter is computed for the entire four- 
dimensional spacctime lattice, I compute and plot E for 
each time slice. Traditional lattice gauge theory is recov- 
ered if a(t) = 1 is used in the equations above, in which 
case E does not have time-dependence. 

For brevity the discussion above is specific to the Z 2 
gauge group, as this is the simplest model to work with 



computationally, but produces the novelties of introduc- 
ing a scale factor. However, the framework for including 
a scale factor in a pure lattice gauge theory discussed 
above works in general for Zn,U(1) and non-Abelian 
gauge groups SU(N). 

The partition function Z still contains a large number 
of terms, 2 640000 for a 20 4 Z 2 lattice, and direct evalua- 
tion on the computer is not feasible. The Monte Carlo 
method is used where a sequence of configurations is gen- 
erated which estimates the full sum, at thermal equi- 
librium at inverse temperature /3, using the heat bath 
method This procedure selects a new value Uij for 
each spin variable in a stochastic manner with proba- 
bility distribution proportional to the Boltzmann factor 
e - ^ , conforming to the condition of detailed balance, 
assuring that the procedure transforms an ensemble in 
equilibrium into itself. 

After setting up the lattice gauge framework, I dis- 
cuss the cosmological background used in the simula- 
tions. In accordance with astrophysical measurements, 
on cosmological scales our Universe is described by the 
spatially flat, homogeneous and isotropic FLRW metric 
ds 2 = a 2 {t)dY? - dt 2 , where dY? = d 2 x + d 2 y + d 2 z. 

With the FLRW metric the Einstein equations, de- 
scribing the interaction of spacetime and matter, reduce 
to the Friedmann equations, a set of second order differ- 
ential equations which can be solved for a(t) if an equa- 
tion of state p = p(p) is specified, where p is pressure 
and p is energy density. In cosmology the perfect fluid 
equation of state p = wp is used, characterized by the di- 
mensionless number w. Different values of w correspond 
to different source types. 

In this Letter I work with three source types commonly 
considered in cosmology: matter, radiation and lambda. 
Table U gives w and the form of the scale factor a(t) 
arising from the Friedmann equations for these source 
types. 



TABLE I. Scale factor for standard source types. 



source 


w 


a(t) 


matter 





■1 

t~ 


radiation 


1/3 


1 


lambda 


-1 


exp(t) 



The simulations are carried out on a 40 x 20 3 hypercu- 
bic lattice using periodic boundary conditions. I let the 
spacetime expand until 7/, = 20 and then symmetrically 
contract it back by setting a(Th + 1) = a(T/j — t) so that 
it can be wrapped onto itself at T = 40. This gives an 
effective lattice size of 20 4 , as the thermodynamic quan- 
tities of the two halves are mirror images of each other. 
Only the first, expanding half is shown in the following 
figures. 
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I run simulations using different /3s, starting at (3 = 
to /3 = 0.5 at A/3 = 0.025 intervals. Before each /3-run, 
the lattice is initialized to Uj,j = 1. I experimented with 
different initialization strategies, such as completely dis- 
ordered states, but these did not affect the overall results. 
For each /3-run, the entire lattice is updated 100 times us- 
ing the Monte Carlo heat bath algorithm, a value which 
was found to give good convergence results. The size of 
the lattice, /3-resolution and number of Monte Carlo it- 
erations used is only limited by computational resources. 
Comparing my results with past results [l[ shows that 
these values arc sufficient for good convergence and un- 
derstanding the thermodynamic properties of the system. 
With these values, a simulation run takes a few minutes 
on a standard desktop computer. 

Figure [T] shows the results of simulation runs for the 
three source types considered. A pixel on each of these 
2D maps is the average action E of the 4D lattice for a 
simulation run at a given /3 (row), at each time step t 
(column) after the 100 Monte Carlo iterations have com- 
pleted. Since I am working in 4D spacetime, the entire 
row is extracted at once at the end of each /3-run. For 
a fixed /3 below (3 C = 0.44, the Universe starts off in a 
completely disordered state, at some /3 dependent time 
reaches a phase transition and then becomes completely 
ordered. Examining the three cases (a) - (c), the transi- 
tion occurs sooner for cosmologies where the scale factor 
a{t) grows faster, lambda being the fastest and radia- 
tion the slowest. Above /3 C , the Universe is ordered at 
all times and no phase transition occurs inside the space- 
time lattice. The critical value /3 C = \ ln(l + \/2) w 0.44 
is known from the self-duality of the Z2 theory [8| . 

The black line on all three plots corresponds to 
f3 cff (t) = (3 [a 5 (t) + a 7 (t)] /2 = /3 C , and is displayed to 
show that although different scale factors are employed 
in the three simulations, there is a common theme. The 
presence of the scale factor introduces the effective in- 
verse temperature /3 c s(t), whose form is explained by in- 
specting the expression for the per plaquette action S\j. 
The metric factors in S\j give either a 5 (t) or a 7 (t), the 
former for timclike plaqucttcs where cither fi or v is 0, the 
latter for spacelike plaquettes where neither \i or v is 0, 
and there are an equal number of timelikc and spacelike 
plaquettes in a hypercubic lattice. 

To verify that the cosmological model with a scale 
factor is statistically equivalent to one with a time- 
dependent effective /3, I have performed simulations on 
the same lattice, removing the scale factor but letting 
/3 = /3 ff(i) vary with time according to the functional 
forms given in Table HI In this approximation both time- 
like and spacelike per-plaquette actions are computed at 
the same f3 e s(t) inverse temperature at each time slice. I 
have subtracted these maps from the cosmological maps 
to examine the error of this estimation. The results are 
shown on Fig. [5] and demonstrate good agreement, with 
the largest deviation occurring at the phase transition 
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FIG. 2. Difference in the order parameter E for an effec- 
tive 20 4 lattice between Monte Carlo simulations where a 
cosmological scale factors was used versus a time-dependent 
f3 e ff(t) = P [a 5 (t) + a 7 (t)] /2. The black line corresponds to 
P e s{t) = Pa- The plots demonstrate that the statistics of the 
system is well approximated by the time-dependent effective 
P. Since P = 2 /go, this corresponds to a time-dependent 
coupling constant in the gauge theory. 



around the /3 c ff(t) = /3 C line, when the system in under- 
going the largest fluctuations. 

To further examine the /3 c ff(t) approximation, I have 
run high resolution simulations at fixed /3 = 0.25 for the 
radiation case, with both the cosmological scale factor 
and the effective inverse temperate on a lattice with ef- 
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FIG. 3. E for high resolution Monte Carlo runs for cosmology 
(crosses) and time-dependent effective f3 (squares), for fixed 
/3 = 0.25. In both cases the scale factor a(t) for radiation is 
used. The vertical line corresponds to t — 6.58 where j3 e s{i) = 

fective size 200 x 20 3 using 100 Monte Carlo iterations. 
Results are shown in Fig. [3] and show excellent agree- 
ment between the two models. At t = the order pa- 
rameter starts out at E ps 0.75, the value for E when 
performing a traditional Zi 4D simulation at f3 = 0.25. 
The vertical line is the critical j3 c value at which the phase 
transition is predicted to occur. As first shown by Creutz 
Q j Fig- G2 demonstrates that the Zi model is first-order 
discontinuous with its sharp drop at j3 c . 

The traditional Z 2 model shows hysteresis in E when 
a thermal cycle is simulated on the system, which cor- 
responds to first increasing /3 from to above /3 C , and 
then back to 0. For the system studied in this letter, 



the effective (5 is a monotonic function of a(t) in the rel- 
evant region a(t) > 1. The simulation is set up in such 
a way that a(t) increases in the first half of the lattice, 
for t < Th, and then decreases symmetrically for t > Th, 
so periodic boundary conditions can be used. However, 
comparing the two sides of the lattice I found no hys- 
teresis effect when flipping over E about Th- Hysteresis 
occurs in the original model because as (3 is increased (or 
decreased) for each iteration, the lattice is in the state 
of the previous iteration, and this dependence results in 
different curves for the two /3 directions. A quick modi- 
fication to the cosmological simulation verifies that hys- 
teresis can be reintroduced by setting the [Zy variables 
at time slice ij+i equal to the corresponding variables at 
ti when performing the simulation. 

I plan to study this modified system, along with the 
more general Zjv, U(l) and SU(N) gauge groups in the 
future. 
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